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Abstract 



The infinite-range-interaction Ising spin glass is considered in the presence 
of an external random magnetic field following a trimodal (three-peak) dis- 
tribution. Such a distribution corresponds to a bimodal added to a proba- 
bility pq for a field dilution, in such a way that at each site the field hi obeys 
P{hi) = p + 5(hi — ho) J t-po5(hi)+p-5(hi + ho). The model is studied through 
the replica method and phase diagrams are obtained within the replica- 
symmetry approximation. It is shown that the border of the ferromagnetic 
phase may present, for conveniently chosen values of po and ho, first-order 
phase transitions, as well as tricritical points at finite temperatures. Anal- 
ogous to what happens for the Ising ferromagnet under a trimodal random 
field, it is verified that the first-order phase transitions are directly related 
to the dilution in the fields: the extensions of these transitions are reduced 
for increasing values of po- Whenever the delta function at the origin be- 
comes comparable to those at hi = ±/to, first-order phase transitions disap- 
pear; in fact, the threshold value Pq, above which all phase transitions are 
continuous, is calculated analytically as p^ = 2(e 3 / 2 + 2) _1 m 0.30856. The 
ferromagnetic boundary at zero temperature also exhibits an interesting be- 
havior: for < po < Pq, a single tricritical point occurs, whereas if po > Po 
the critical frontier is completely continuous; however, for po = Pq, a fourth- 
order critical point appears. The stability analysis of the replica-symmetric 
solution is performed and the regions of validity of such a solution are iden- 
tified; in particular, the Almeida-Thouless line in the plane field versus 
temperature is shown to depend on the weight po- 
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1. Introduction 



Among disordered magnets [y, spin glasses [@-|4J] and ferromagnets in the presence 
of random fields may be singled out as two of the most puzzling and controversial 
systems in condensed matter physics. 

The random- field Ising model (RFIM), introduced by Imry and Ma ||, has concen- 
trated a lot of interest after the identification of its physical realizations. Probably the 
most important physical conception of the RFIM comes out to be a diluted Ising antifer- 
romagnet in the presence of a uniform magnetic field [10|,p|. Since then, many diluted 



antiferromagnets have been investigated, in such a way that systems like Fe x Zni_ x F2 and 
Fe x Mgi_ x Cl2 are nowadays considered as standard experimental realizations of the RFIM 



12| . [T3|] . From the theoretical point of view, many important ingredients remain unknown. 
At the mean-field level, it is well known that different probability distributions for the 
random fields may lead to distinct phase diagrams, e.g., a Gaussian probability distribu- 



tion yields a continuous ferromagnetic-paramagnetic boundary [rjj, whereas for a bimodal 



distribution, this boundary exhibits a continuous piece at high temperatures ending up at 
a tricritical point, which is followed by a first-order phase transition at low temperatures 



[ 375 1 . Such a contrast in the mean-field phase diagrams of the RFIM with the bimodal 



and Gaussian probability distributions has been proven rigorously ||16 . Indeed, Aharony 



1 15] argued that whenever an analytic symmetric distribution for the fields presents a 
minimum at zero field, one should expect a tricritical point and a first-order transition for 
sufficiently low temperatures. Further studies of the RFIM at the mean-field level have 



considered a trimodal (three-peak) distribution (pTTLlIB 



P(hi)=p + 6(hi-ho)+po5(hi)+p-5(hi + ho) , (1.1) 
in its symmetrical form, i.e., p + = p~ = |(1 — po). Such a distribution, which may 



be interpreted as a bimodal added to a dilution in the fields with probability po 13 > ^ s 
expected to mimic better real systems than its bimodal counterpart. It was shown that 
the field dilution plays an important role in what concerns the presence of the tricritical 
point: distinct analyses lead to slightly different estimates for the threshold value, above 
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which the tricritical point disappears (whereas the analysis of Mattis [|17] shows that the 



tricritical point vanishes for po > 0.25, according to Kaufman et al. [18] such a behavior 
should occur for po > 0.24). Whether the features in the mean-field phase diagrams of 
the RFIM should prevail on short-range-interaction models, represents a point which has 
attracted a lot of interest ]|l9| -|2"2"|j . For the three-dimensional RFIM, recent Monte Carlo 
simulations detect a jump in the magnetization but no latent heat, for both bimodal [13] 



and Gaussian [20] distributions, whereas high-temperature series expansions [21] and a 
zero-temperature scaling analysis |22| find a continuous transition for both distributions. 
However, in four dimensions the same zero-temperature analysis |22] leads to a first- 
order transition in the bimodal case and a continuous one for a Gaussian distribution, in 
agreement with the mean- field predictions. Apart from that, the low-temperature phase 
of the RFIM, in finite dimensions, may present a nontrivial structure, with a complicated 
free-energy landscape, as suggested by perturbative analyses ]p3|p4 . 



The Ising spin-glass (ISG) problem became, nowadays, one of the most controversial 
issues in the physics of disordered magnets. Its mean-field theory, based on the solution 
of the infinite-range-interaction model, the so-called Sherrington-Kirkpatrick (SK) model 
||25|| , presents a quite nontrivial behavior. The correct low-temperature solution, as pro- 
posed by Parisi [2~B|, consists of a continuous order-parameter function (i.e., an infinite 
number of order parameters) associated with many low-energy states, a procedure which 
is usually denominated as replica-symmetry breaking (RSB). Furthermore, a transition 
in the presence of an external magnetic field, known as the Almeida-Thouless (AT) line 



27 1 , is found in the solution of the SK model: such a line separates a low-temperature re- 



gion, characterized by RSB, from a high-temperature one, where a simple one-parameter 
solution, denominated as replica-symmetric (RS) solution, is stable. The validity of the 
results of the SK model for the description of real (short-range-interaction) systems rep- 
resents a very polemic question pi. The rival theory is the droplet model | 2~B| , based 
on domain- wall renormalization-group arguments for spin glasses [ P§PU| . According to 
the droplet model, the low-temperature phase of any finite-dimensional short-range spin 
glass should be described in terms of a single thermodynamic state (together, of course, 
with its corresponding time-reversed counterpart), i.e., essentially a RS-type of solution. 
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Obviously, the droplet model becomes questionable for increasing dimensionalities, where 
one expects the existence of a finite upper critical dimension - believed to be six for 



the ISG jST] - above which the mean-field picture should prevail. Recent analyses of 
short-range ISG on diamond hierarchical lattices (on which the Migdal-Kadanoff renor- 
malization group is exact) has found evidences of the droplet picture [3_^]; however, the 
applicability of such lattices for the description of ISG on Bravais lattices is doubtful 
|33|,34]. Numerical simulations are very hard to be carried for short-range ISG on a cubic 



lattice, due to large thermalization times ||33|| ; as a consequence, no conclusive results in 
three-dimensional systems are available. However, in four dimensions the critical tem- 
perature is much higher, making thermalization easier; in this case, many works claim to 
have observed some mean-field features [55]. 

From the theoretical point of view these two problems (RFIM and ISG), have been, in 
most of the cases, studied in separate, with a few exceptions [SB-Eo| . However, many di- 



luted antiferromagnets, like Fe x Zni_ x F2 [il[] and Fe x Mgi_ x Cl2 |42| , f43|| , are able to exhibit, 
within certain concentration ranges, random-field, spin-glass or both behaviors. For the 
Fe x Zn!_ x F 2 , one gets a RFIM for x > 0.40, an ISG for x < 0.24, whereas for intermedi- 
ate concentrations (0.24 < x < 0.40) one may observe both behaviors depending on the 
magnitude of the applied external magnetic field [RFIM (ISG) for small (large) magnetic 



fields], with a crossover between them; this latter effect was observed in Feo.31Zno.69F2 [41]. 
Certainly, such properties are expected to be properly explained only if one considers a 
model which takes into account both spin-glass and random-field ingredients. Indeed, the 
crossover observed in Fe .3iZn .69F 2 was also found in the study of the SK model under a 
Gaussian random field ||38|| . On the other hand, the study of the SK model in the presence 
of a bimodal random field produced interesting results, with first-order phase transitions 
and tricritical points p£]; such results may be relevant for explaining the first-order phase 
transitions observed in FexMgx.xCL? [ p~3|| . 

In the present work we study the SK model in the presence of a random field following a 
trimodal probability distribution [see Eq. (1.1)]. In addition to that, one may interpolate 
between the bimodal distribution and a behavior which is qualitatively analogous to the 
Gaussian one, since by monitoring the delta function at the origin, one is able to control 
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the presence of tricritical points. In the next section we define the model and, through the 
use of the replica method, we find its free-energy density, equations of state and equations 
for the validity of the RS solution. In section 3 we exhibit and discuss the phase diagrams 
of the model. Finally, in section 4 we present our conclusions. 



2. The Model and Replica Formalism 



The mean-field theory of the ISG is usually formulated as a set of N spins, each of 
them interacting with all others [a total of ^N(N — 1) interactions], known as the SK 



model [25]. The SK model in the presence of an external random magnetic field may be 
defined in terms of the Hamiltonian [B8.B9], 



n = - y £j ij s i s j - y £h i s i , (2.i) 

(ij) 1 

where Si = ±1 , with i = 1,2, ■ • • , iV, and the interactions are infinite-range-like, i.e., 
the sum J2u,j) applies to all distinct pairs of spins. The coupling constants {Jij} and the 
random fields {hi} are quenched variables, following independent probability distributions, 

/ N \^ 
P(J - )= v2^J 6XP 

with P(hi) given by Eq. (1.1) (p+ + po + P- = 1)- Let us, for the moment, keep the 
trimodal probability distribution in its general form of Eq. (1.1); later on, we will see that 
the ferromagnetic boundary does not exist for p + p~, and so, in such a case, we will be 
restricted to the symmetrical form p + = p_ = |(1 — po). It should be mentioned that the 
above randomnesses ({Jij} and {hi}) are usually correlated in real systems; herein for the 
sake of simplicity, we shall consider two independent probability distributions. Therefore, 
for a given realization of bonds and site-fields, ({Jij}, {hi}), one has a corresponding free 
energy, F({Jij}, {hi}), such that the average over the disorder, [ jj^, may be performed 
as independent integrals, 
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N 
2J2 



j. — 11 

13 N 



•Jn 



(2.2) 



MUvh {hi})h,h = J WJijPiJij)} J{[dhiP{hi)\ FttJtj}, {h t }) . (2.3) 



(U) 



The usual procedure consists in applying the replica method [^Q, in such a way as to 
let the free energy per spin as, 

-fif= lim 1 [In Z({J tl },{h t })] JA 

N^oo iV 

= hm lim -La^U-l) , (2.4) 

AT— >oo n^O iVTi 



where Z n is the partition function of n copies of the system defined in Eq. (2.1) and 
(3 = l/T (we work in units ks = 1)- Standard calculations lead to 

(ST) 2 1 

13 f = -StL_L + Jim _ min g ( m <*^ ; (2.5) 
4 n->o n 



where 



p ln Tr a exp(7Y° // ) - p_ln Tr a exp(H e// ) , (2.6a) 



^/z = /^o E™° ^ + J ) 2 E ^ ^ ± Pho E ^ > (2-6b) 

a (a/3) a 

= PJo E m " 5 " + (0 J ) 2 E ^ sast3 ■ ( 2 - 6c ) 

a (a/3) 



In the equations above, the sum indexes a and (3 (a, (3 = 1,2, ■ ■ ■ ,n) are replica labels 
and J2(a/3) denote sums over distinct pairs of replicas. 

The extrema of the functional g(m a ,q a ^) give us the equilibrium equations for the 
magnetization and spin-glass order parameters, respectively, 



m a =p + (S a ) + + p (S a ) +P-(S a ). , (2.7a) 
q^=p + (S a S' 3 ) + +p Q (S a S%+p_(S a S^^ , (2.7b) 

where ( )± and ( )o refer to thermal averages with respect to the "effective Hamiltonians" 
Ti-fff and Ti-^ff * n Eqs. (2.6b) and (2.6c), respectively. 
If one assumes the replica-symmetry (RS) ansatz p5| . 



m a = m , Va ; q a(3 = q , V(a/3) , (2.8) 
the free energy per spin (Eq. (2.5)) and the equilibrium conditions (Eqs. (2.7)) become 

/?/ = -^(1 - g) 2 + ^ m 2 - p + J P,ln(2cosh£+) 



where 



and 



Po 



J Dzln(2cosh£°) ~ P- J £^ln(2cosh£-) , (2.9) 



m= p+ J X>,ztanh£ + + Po J Dztanh£° + p^ J £>ztanh£" , (2.10) 
q= p + J £>ztanh 2 £+ + Po J ^tanh 2 «^°+ p^ J £>ztanh 2 «T , (2.11) 



oo 



2tt 



Vz--- = r {^-Y dz exp(-^ 2 /2) •■■ , (2.12) 



£± = (3J m + (3Jq 1/2 z±f3h , (2.13a) 
£° = (3J m + (3Jq 1/2 z . (2. 13b) 



Although the spin-glass order parameter (Eq. (2.11)) is always induced by a nonzero 
random field (po < 1), it may still contribute to a nontrivial behavior; this is provided by 
the instability of the RS solution. Such an instability occurs at the AT line p7] . 
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= p + Vz sech 4 £ + + p Vz sech 4 £° + p_ Vz sech 4 ^" , (2.14) 



which may be obtained through the simultaneous solution of Eqs. (2.14), (2.10) and 
(2.11). 

In the next section we shall consider the phase diagrams of the model and the regions 
of instability of the RS solution, worked out from Eqs. (2.9)-(2.14). 



3. Results and Discussion 

Let us first consider the case Jq = 0; one may easily see that the only nontrivial 
behavior in this case is given by the AT instability in the plane magnetic field versus 
temperature, which may now be obtained from the solution of Eqs. (2.11) and (2.14). 
The integrals involving £~ may be easily transformed through the change of variables 
z — > —z , in such a way that the AT line may be obtained by solving the set of equations, 



J 



1 -p ) J Vz sech 4 (/5Jg 1/2 z + (3h ) + p J Vz sech\[3Jq 1/2 z) , (3.1a) 



q= (1-po) J Vz tanh 2 (J3 J q 1/2 z + (3h ) + p J Vz ttmh 2 ([3Jq 1/2 z) . (3.1b) 

It should be pointed out that the equations above are valid for arbitrary values of the 
weights in the probability distribution of Eq. (1.1), with p + + p- = 1 — po ; although the 
AT line changes with field dilution, it is no altered under a field inversion. The AT lines 
in the plane magnetic field versus temperature are exhibited in Fig. 1, for typical values 
of po- Clearly, the AT line for the bimodal distribution (p = 0) |39] is identical to the 
one of the SK model in the presence of a uniform magnetic field |27| , due to the property 
of invariance under field inversion. For < p < 1, one may calculate analytically the 
behavior of the AT line in the low-field regime (T = J), 



T 
1 



3(1 -Po) 
4 



1/3 / s2/3 

7 



(3.2) 



which leads to a slightly modified amplitude, but the same low-field exponent of the 
standard AT line If one considers po ~ 0, the low-temperature behavior of the AT 
line may be easily calculated, 



T 4 1 



J 3J2^ 



1 -p )exp - 



2 J 2 



(3.3) 



which exhibits the usual exponential decay ]27], but with a shift towards higher temper- 
atures for increasing values of po- In all other situations, the AT lines were calculated 
by solving numerically Eqs. (3.1). One notices that for high values of po, the integrals 
multiplying po in Eqs. (3.1) contribute significantly, in such a way that the AT lines 
become slightly independent of ho, for ho large enough, as shown in Fig. 1. 

From now on, we will be restricted to Jo > 0; in this CclSG, clS far as RS is concerned, 
if p + 7^ p- Eqs. (2.10) and (2.11) yield nonzero magnetization and spin-glass order 
parameters, leading to trivial behavior. Therefore, for the rest of this paper we will 
concentrate on a symmetrical trimodal distribution, i.e., p+ = P- = |(1 — Po)- In this 
case, the random field still induces the parameter q, leading to no spontaneous spin-glass 
order (like the one found for the SK model in the absence of external field [p5]). Therefore, 
the only possible phase transition within the RS approximation is the one associated with 
the magnetization, similarly to what happened in the case of the bimodal distribution 
|39| . Hence, two phases are possible, namely, the ferromagnetic (m/0 , q ^ 0) and the 
independent (m = , q ^ 0) ones. Although in the RFIM this latter phase is usually 
denominated of paramagnetic, in the present problem, within the RS approximation, we 
shall keep the nomenclature independent, for reasons which will become clear soon. 

The critical frontier separating these two phases may be found by solving the equi- 
librium equations, (2.10) and (2.11); in the case of first-order phase transitions, we shall 
make use of the free-energy per spin [Eq. (2.9)] as well. Expanding Eq. (2.10) in powers 
of m one gets, 

m = Ai(q) m + A 3 (q) m 3 + A 5 (q) m 5 + 0(m 7 ) , (3.4) 
where the coefficients depend on q [which on its turn, depends on m through Eq. (2.11)]. 
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Expanding Eq. (2.11) in powers of to, 



(pjo) r 2 . 4 . 
g = go i _ ijjf r m ^ m ^ ' ^ ^ 

with 

Po)(l - 4p+ + 3p+) +p (l - 4p? + 3p°) , (3.6) 

J Vz ttmh 2k ({3Jql /2 z + (3h ) , (3.7a) 
J Vz tanh 2fe {(3 J q^z) , (3.7b) 

where q is independent of to, corresponding to the solution of Eq. (2.11) with to = 0. 
Substituting the above results into Eq. (3.4), one gets the m-independent coefficients of 
the power expansion, 

A[ = pj [l - (1 - p )p+ - pop ,] , (3.8a) 



r = (i- 

pt = 
pI = 



i + 2(/?j) 2 r" 
i - (pjy r _ 



(3.8b) 



A' = 



-7 



( W 

30 



8((3J) 2 T + 36([3J) 4 T 2 + lh((3Jf T 3 
1 - {(3 J) 2 r 



(3.8c) 



where 



7 = (1 - po)(-4 + 34p+ - 60p+ + 30p+) + p (-4 + 34p? - 60p° + 30p° 



(3.9) 



The critical frontier may be determined using standard procedures, as described below. 

(i) For continuous phase transitions, A[ = 1 and A' 3 < . 

(ii) A first-order phase transition occurs whenever A\ = 1 and A' 3 > 0; the proper 
critical frontier should be found, in this case, through a Maxwell construction, i.e., by 
equating the free energies of the two phases. 
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(iii) When both types of phase transitions are present, the continuous- and first-order 



critical frontiers meet at a tricritical point [45], which defines the limit of validity of the 
series expansions; beyond the tricritical point the magnetization is discontinuous. The 
location of such point is determined by setting A[ = A' 3 = 0, with the condition A' 5 < 
satisfied. 

In Figs. 2-4 we show three qualitatively distinct ferromagnetic boundaries of the 
present problem, for a typical value of po (po = 0.3), compared with those of the bimodal 
probability distribution (p = 0). In Fig. 2 there is a single point along the ferromagnetic 
boundary at which A' 3 = 0; such a point may not be considered as tricritical, since there 
is no first-order phase transition. However, for any value of ho greater than those of Fig. 
2 [h /J = 0.9573 (po = 0) and ho/ J = 1.53526 (po = 0.3)], one gets first-order phase 
transitions, and at least one tricritical point. In Fig. 3 we show situations where two 
tricritical points appear along the ferromagnetic boundary; we have verified that, for a 
fixed value of po, such a behavior occurs within a narrow interval of ho- In Fig. 4 a 
single tricritical point emerges, separating a continuous boundary (high temperatures) 
from a first-order critical frontier (low temperatures). From such phase diagrams, one 
notices that the main effect of the field dilution is to push the tricritical points towards 
lower temperatures, i.e., the temperature range over which the first-order transitions occur 
decreases. 

As mentioned before, although the spin-glass order parameter is always induced by 
the random field, it may still exhibit interesting behavior, associated with the instability 
of the RS solution. The AT instabilities, given by the solution of Eqs. (2.10), (2.11) 
and (2.14) with p + = p_ = |(1 — p ), yields two distinct lines in the phase diagrams of 
Figs. 2-4, depending on whether one is inside the independent phase (m = 0), or in the 
ferromagnetic (m ^ 0) one. In the former case, the AT line is a straight line (independent 
of Jo), whereas in the latter, it presents the usual decrease with temperature for increasing 
values of Jo, in such a way that for low temperatures one gets the exponential decays, 



T 
1 



4 1 



1 



3J2n 2 



1 - po) exp 



(Jo + hp) 
2J 2 



+ p exp 



J 



2 J 2 



12 



+ 2^ -Po) ex P 

Herein we shall adopt the usual criteria for the identification of the regions where RS is 
stable and those throughout which a RSB procedure is necessary |3||§]. The two regions 
with zero magnetization will be associated with the paramagnetic (high temperatures) 
and spin-glass (low temperatures) phases, whereas those with nonzero magnetization will 
be associated with the ferromagnetic (high temperatures) and mixed-ferromagnetic (low 
temperatures). The several phases exhibited in our phase diagrams are identified as: 



(Jo ~ hp) 
2 J 2 



(3.10) 



Paramagnetic (P) 


(m 







q : 


RS) ; 


Spin-Glass (SG) 


(m 







q 


RSB) ; 


Ferromagnetic (F) 


[m 







q 


RS) ; 


Mixed Ferromagnetic (F') 


[m 







q 


RSB) . 



It should be mentioned that the present low-temperature results are questionable inside 
the phases F' and SG, due to the instability of the RS solution; in particular the point 
for po — 0.3 where A' 3 = in Fig. 2, as well as the low-temperature tricritical points of 
Fig. 3 may completely disappear under a RSB procedure. However the high-temperature 
tricritical points, like those of Figs. 3 and 4, are inside the region of stability of the RS 
solution and will persist under more general treatments; we believe that such points are 
reminiscent of the tricritical point of the bimodal RFIM. 

The two AT lines mentioned above usually meet at a continuous ferromagnetic bound- 
ary; however, these lines do not match each other across first-order phase transitions 
39| , |40|^6| : there is a small (but finite) gap between them in Figs. 3 and 4. 



Let us now investigate the ferromagnetic boundary at zero temperature; for T = 
the spin-glass order parameter is trivial (q = 1), in such a way that one gets for the free 
energy and magnetization, 



J 2 h 
J = 7T m ~ TTv "Po) 



fj m + h \ (J Q m-h Q 

erf -= — — erf = — 

\ JV2 \ Jy/2 . 
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J 



l-po) lexj) 



(J Q m + h Q f 



2J 2 



exp 



( J m - h c 
2J 2 



2J 



Po \ exp 



(Jomf 
2J 2 



m 



erf 



' J m + /i c 
4 Jv^ 



+ erf 



J m — ho 



/ J m N 

+Poerf (vs, 



(3.11a) 



(3.11b) 



Using a similar procedure as the one for finite temperatures, one may expand Eq. (3.11b), 

m = ai m + a 3 m 3 + a 5 m 5 + 0(m ) , (3-12) 

where, 



ai 



0.3 



a 5 



2 Jo 

7T J 



(1 - p ) exp 



''0 



2 J 2 



+ Po 



1 

6 



V 3 



(1-Po) [ f 2 -l exp 



/r 



2 J 2 



Po 



120 



7 \ 5 

J 



;i _ Po) (M_ 6 | + 3 )exp(-§- 



3pc 



(3.13a) 



(3.13b) 



(3.13c) 



The critical frontier separating the phases F' and SG is shown in Fig. 5 for typical values 
of po- One notices that the effect of the weight po is to favour the continuous line, along 
which ai = 1 with < 0, i.e., 



Jo 



J V2 Po + (l- Po )exp(-/ i 2 /2J 2 ) 



(3.14) 



while decreasing the extension of the first-order transition line. For small values of po these 
two lines meet at a tricritical point, obtained by solving the equations a± — 1, 03 = 0, 
with the condition 05 < 0; within the analysis for finite temperatures, this corresponds 
to the situation where the lower-temperature tricritical point (cf. Fig. 3) hits the zero- 



temperature axis. If po = such an effect occurs at [39 



14 



hn Jo fife 

-1 = 1 ; ^=^2.0664 



(3.15) 



We verified that for < p < Pq (where p^ will be defined below), such a set of equations 
presents two solutions, although only one of them represents a tricritical point, satisfying 
a 5 < 0. By increasing p Q inside this range, we noticed that such solutions get closer and 
colapse for p = p* . We calculated analytically p* = 2(e 3 / 2 + 2) _1 0.30856, at which a 
fourth-order critical point j|7| (characterized by a\ = a% = = 0, with a-j < 0) occurs at 

^ = V3w 1.73207 ; 4 = — (e 3/2 + 2) w 2.70786 . (3.16) 
J J 6 

The value represents a threshold of po, above which there are no first-order transitions 
for any temperature T > 0. For p > Pq the second-order critical frontier of Fig. 5 
approaches an asymptote for large values of ho] indeed, when p — > 1 the zero-temperature 
ferromagnetic boundary approaches a straight line at Jo/ J = yV/2 [see Eq. (3.14)], 
characteristic of the SK model in zero field [pB]| . 

It should be mentioned that the finite-temperature vestigial points where A' 3 = 0, like 
the ones in Fig. 2, are qualitatively different from the fourth-order critical point found 
for po = p^ at zero temperature, even though both situations represent thresholds for the 
occurrence of tricritical points. In the former case, A' 5 < 0, whereas in the latter, A' 5 = 0. 
In Fig. 6 we exhibit the behavior of the coefficients A' 3 and A' 5 , for temperatures along 
the ferromagnetic frontier, for the case (b) of Fig. 2, i.e., p = 0.3 (h / J = 1.53526), and 
Po = Pq (ho/ J = One clearly sees that the fourth-order critical point only shows up 
at zero temperature; its parameters, as defined in Eq. (3.16), correspond to the situation 
where the vestigial point of Fig. 2 collapses with the zero-temperature axis. 

If < po < Pq, it is always possible to obtain first-order phase transitions by conve- 
niently choosing the value of h . In Fig. 7 we exhibit the ranges of po and h /J throughout 
which first-order phase transitions and tricritical points are possible along the ferromag- 
netic boundary. In region (a), first-order phase transitions are conceivable at finite and 
zero temperature, with a single tricritical point (at finite temperatures): typical examples 
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are shown in Fig. 4. Throughout a very narrow range [region (b)] two tricritical points 
appear and the first-order phase transition occurs only for finite temperatures: typical 
examples are exhibited in Fig. 3. The region (b) is delimited by characteristic values of 
(po, ho/ J): (i) the threshold for ho/ J smaller corresponds to the set of points satisfying 
A' 3 = 0, but with no first-order phase transition (e.g., the vestigial points shown in Fig. 
2); (ii) the delimiter for ho/ J larger corresponds to the coordinates of the tricritical points 
at zero temperature. The vertical line in Fig. 7 is for po = Po, defining [together with 
the delimiter (i) of region (b)], the range throughout which the ferromagnetic boundary 
is always continuous [region (c)]. 



4. Conclusion 

We have studied the Sherrington-Kirkpatrick spin glass in the presence of random 
fields {hi}, following a trimodal (three-peak) probability distribution, which corresponds 
to a bimodal plus a probability p for field dilution, i.e., P(hi) = p + 8(hi — ho) + PoS(hi) + 
P-5{hi + ho). We have used the replica method and the phase diagrams were obtained 
within the replica-symmetry approximation. The boundary of the ferromagnetic phase 
exhibited an interesting behavior, with the presence of first-order phase transitions and 
tricritical points: within certain ranges forpo and ho, a single or two tricritical points were 
encountered. We have shown that the first-order phase transitions are directly affected 
by the dilution in the fields, in such a away that the extension of such lines are reduced 
by increasing p . In fact, there is a threshold value, p^ = 2(e 3 / 2 + 2) _1 « 0.30856, above 
which the ferromagnetic boundary is always continuous. Such effects may be reminiscent 
of those occurring within the mean-field theory of the Ising ferromagnet in the presence of 
trimodal random fields: the single tricritical point that appears in the case of a bimodal 
distribution JT51 is washed way by the presence of the delta at the origin, whenever po 



becomes greater than a certain value [17,18 



At zero temperature, if < po < Po, the ferromagnetic critical frontier exhibits a single 
tricritical point, with a first-order phase transition at high values of h . By increasing p , 
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the first-order line gets reduced and, for p = p^, a fourth-order critical point is observed; 
for p > Po, the ferromagnetic boundary is always continuous. 

Although the spin- glass order parameter is induced by the random field (po < 1), it 
may still contribute to a nontrivial behavior, in what concerns the stability of the replica- 
symmetric solution. We have calculated the regions of instability of such a solution, 
leading to the identification of two low-temperature phases, namely, the spin-glass and 
mixed ferromagnetic ones. Besides that, the Almeida-Thouless line in the plane field 
versus temperature was shown to depend on the weight po, with different amplitudes (but 
the same exponent) in the low-field regime, and qualitatively distinct high-field behaviors. 

We have verified that whenever the ferromagnetic boundary presents both continuous 
and first-order transition lines meeting at a single finite-temperature tricritical point, 
such a point is located inside the region of stability of the replica-symmetric solution, 
and it will not be removed by a replica-symmetry-breaking procedure. However, when 
two tricritical points occur along the ferromagnetic boundary, at least one of them (the 
one at low temperatures) appears inside the unstable region, and its existence may be an 
artifact of the replica-symmetric solution. 

The applicability of the present results in the description of real systems obviously 
depends on the survival of the mean-field characteristics in the respective short-range- 
interaction versions of Ising spin glasses and the Ising ferromagnet in the presence of a 
random field. However, the trimodal distribution employed herein is expected to mimic 
better real systems than the bimodal distribution itself. Although we are not aware 
of experimental observations that match with our results, we believe that the diluted 
antiferromagnet Fe x Mgi_ x Cl2 is a good candidate, since, for conveniently chosen dilutions, 
it may exhibit first-order phase transitions [Tj| , as well as a crossover from first- to second- 



order behavior [44 
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Figure Captions 



Fig. 1: The AT lines, for the SK model in the presence of a trimodal random field, in 
the plane h versus T (in units of J), for typical values of p . 

Fig. 2: Phase diagram T versus J (in units of J) of the SK model in the presence of 
a trimodal random field with p — 0.3, compared with one of the bimodal case (po = 0), 
for conveniently chosen values of ho- (a) h /J = 0.9573 (po = 0); (b) ho/ J = 1.53526 
(po = 0.3). The ferromagnetic boundaries are continuous, except for the points where 
A' 3 = [cf. Eq. (3.7b)], represented by black squares. These choices signal lower 
bounds for ho, above which first-order phase transitions occur. The phase nomenclature 
is specified in the text, with the low-temperature phases SG and F' delimited by AT 
lines. 

Fig. 3: Phase diagram T versus J (in units of J) of the SK model in the presence of 
a trimodal random field with p = 0.3, compared with one of the bimodal case (po = 0), 
for conveniently chosen values of ho, in such a way as to obtain two tricritical points 
(black circles) along the ferromagnetic boundary, (a) h / J = 0.97 (p = 0); (b) h / J = 
1.558 (po = 0.3). The dashed lines stand for first-order phase transitions. The phase 
nomenclature is the same as in Fig. 2. 

Fig. 4: Phase diagram T versus J (in units of J) of the SK model in the presence of 
a trimodal random field with p = 0.3, compared with one of the bimodal case (po = 0), for 
conveniently chosen values of ho, in such a way as to obtain a single tricritical point (black 
circle) along the ferromagnetic boundary, (a) ho/ J = 1.02 (p = 0); (b) ho/ J = 1.58 
(po = 0.3). The phase nomenclature and line representations are as in Figs. 2 and 3. 



22 



Fig. 5: The zero-temperature phase diagram h versus J (in units of J) of the SK 
model in the presence of a trimodal random field, for typical values of po- If < po < Pq 
one always gets tricritical points (black circles), followed by first-order phase transitions 
for high values of ho. When po — Po, one gets a fourth-order critical point (represented 
by a star). Above the threshold value p^ = 2(e 3 / 2 + 2) _1 0.30856, the critical frontier 
separating the phases SG and F' is continuous. 

Fig. 6: The ordinate represents either the coefficient A' 3 or A' 5 [Eqs. (3.7b) and 
(3.7c), respectively] along the ferromagnetic boundary, for p = 0.3 (h /J = 1.53526) 
(dot-dahed lines) and p = p*$ (ho/ J = (full lines), as a function of temperature. In 
the former case, A' 3 = at T/J = 0.25 (with A' 5 < 0), whereas in the latter, A' 3 = A' 5 = 
at T = 0. 

Fig. 7: Ranges of p and h / J associated with distinct behaviors for the ferromag- 
netic boundary, (a) First-order phase transitions at finite and zero temperatures, with a 
single tricritical point at finite temperatures; (b) Two tricritical points with a first-order 
phase transition for finite temperatures; (c) Continuous phase transitions. 
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